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Abstract 

This paper is concerned with uniqueness for reconstructing a periodic inhomogeneous 
medium covered on a perfectly conducting plate. We deal with the problem in the frame of 
^ I time-harmonic Maxwell systems without TE or TM polarization. An orthogonal relation 

for two refractive index is obtained, and then inspired by Kirsch's idea, the refractive index 
can be identified by utilizing the eigenvalues and eigenfunctions to a quasi-periodic Sturm- 
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■ Liouville eigenvalue problem. 
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■ 1 Introduction 

■ Scattering theory in periodic structures has many applications in micro-optics, radar imaging and 
nondestructive testing. We refer to |22| for historical remarks and details of these applications. 
Consider a time-harmonic electromagnetic plane wave incident on a bi-periodic layer sitting on a 
perfectly conducting plate in M.^. We assume that the medium inside the layer consists of some 
inhomogeneous isotropic conducting or dielectric material, whereas the medium above the layer 
consists of some homogeneous dielectric material. Suppose the magnetic permeability is a fixed 
positive constant throughout the whole space. The material properties of the media are then 
characterized completely by an index of refraction in the layer and a positive constant above the 
layer. The direct scattering problem is, given the incident field and the bi-periodic refractive 
index, to study the electromagnetic distributions, whereas the inverse scattering problem is to 
determine the refractive index from the knowledge of the incident waves and their corresponding 
measured scattered fields. 
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Adopting the Cartesian axis 0x1X2X3 with the xs-axis vertically upwards, perpendicular to 
the plate. If the refractive index is invariant in the X2 direction, the direct and inverse prob- 
lems as indicated above can be dealt with in the TE polarization case where the electric field 
E(x) is transversal to the (xi, 3;3)-plane by assuming E = (0, n(xi, 2:3), 0), or in the TM po- 
larization case where the magnetic field H{x) is transversal to the (xi, X3)-plane by assuming 
H = (0, u{xi,X3),0). In the case of TE polarization, Kirsch [17] has studied the direct scattering 
problem via the variational method, and for the inverse problem, instead of constructing the 
complex geometrical optical solutions as in the Calderons problem (see |19l I28|). he considered 
a class of eigenfunctions to a special kind of quasi-periodic Sturm-Liouville eigenvalue problem. 
Relying on the asymptotic behavior of those eigenvalues, the uniqueness result for the inverse 
problem can be proved once the orthogonal relation for two different refractive indexes has ob- 
tained. See also |25|I26) for the direct and inverse acoustic scattering by periodic, inhomogeneous, 
penetrable medium in the whole M?. Other uniqueness results for reconstructing the profile of a 
bi-periodic perfectly conducting grating can be seen in [21 [5l |6]. 

In this paper, we are mainly concerned with the uniqueness issue for reconstructing the refrac- 
tive index in the framework of time-harmonic Maxwell equations without TE or TM polarization. 
The uniqueness result for the inverse problem in this paper is most closely related in term of 
result and method of argument to Kirsch on the determination of the refractive index in the 
TE polarization. Inspired by |27) and |15| . we obtain an orthogonality relation for two different 
refractive indexes by using a D-to-N map on an artificial boundary on which the tangential elec- 
tric fields are identical for an integral type of incident electric field. It should be remarked that 
the method for constructing geometry optical solutions in |19| \TE[ I27| for non-periodic inverse 
conductivity problems does not work since the solutions are required to be quasi-periodic in the 
periodic case. To reconstruct the refractive index, we follow Kirsch's idea [17] (see also [26j) by 
considering a kind of Sturm-Liouville eigenvalue problems. We shall prove the uniqueness result 
when the index depends only on one direction (xi or X2). However, we expect the result to hold 
in a more general case by constructing special solutions with suitable asymptotic behaviors for 
the Maxwell equations. 

Scattering by bi-periodic structures have been studied by many authors using both integral 
equation methods and variational methods (see, e.g. [T], [4], |12| . |13| . |14| . |16) . |20| and [24]). 
It is known that, for all but possibly a discrete set of frequencies, the direct scattering problem 
has a unique weak solution in the case of bi-periodic inhomogeneous medium in the whole M^, 
of which an absorbing medium always leads to a uniqueness result for any frequency. When the 
refractive index is non-absorbing, uniqueness can be guaranteed in the TE mode if the refractive 
index satisfies an increasing criterion in the X3-direction (|25| IT]). See also [10] and [30] for the 
uniqueness results of more general rough surface scattering by an inhomogeneous medium in a 
half space in the TE or TM mode. In this paper, we assume that the medium inside the layer is 
absorbing so that the uniqueness result for the direct problem holds, implying that the D-to-N 
map T (at the end of Section [3]), which depends on the refractive index, is well-defined. 

The rest of the paper is organized as follows. In the next section we set up the precise 
mathematical framework and introduce some quasi-periodic function spaces needed. In Section 
131 we consider a quasi-periodic boundary value problem (QPBVP) in a periodic cell via the 
variational approach which is used for the study of the inverse problem. Uniqueness and existence 
of solutions to the QPBVP are justified by the classic Hodge decomposition and the Fredholm 
alternative. This leads to the definition of a D-to-N map on an artificial boundary which is 
continuous and depends on the refractive index. In Section HI based on the property of the 
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transparent boundary condition defined on the artificial boundary, we give a solvability result 
of tlie direct scattering problem. In Section [5l we establish a uniqueness result for the inverse 
scattering problem. 



2 Time-harmonic Maxwell equations and quasi-periodic function 
spaces 

2.1 Time-harmonic Maxwell equations 

Let = {{xi, X2, X3) G M'^ I X3 > 0} and assume that is filled with an inhomogeneous, 
isotropic, conducting or dielectric medium of electric permittivity e > 0, magnetic permeabil- 
ity //> and electric conductivity o" > 0. Suppose the medium is non-magnetic, that is, the 
magnetic permeability /x is a fixed constant in M.^ and the field is source free. Then the electro- 
magnetic wave propagation is governed by the time-harmonic Maxwell equations (with the time 
variation of the form e"*"^*, u > 0) 

cuil E -iufiH = 0, cmlH + iuj{e + i-)E = 0, (2.1) 

where E and H are the electric field and magnetic field, respectively. Suppose the inhomogeneous 
medium is 27r-periodic with respect to xi and X2 directions, that is, for all ni,n2 G Z^, 

e{xi + 27rni,X2 + 27rn2,X3) = e(xi,X2,X3), a{xi + 27rni,X2 + 27rn2,X3) = cr(xi, X2, X3). 

Further, assume that e(x) = cq, a = for X3 > 6 (which means that the medium above the 
layer is lossless) and that the inhomogeneous medium has a perfectly conducting boundary 
Tq := {x3 = 0}. Consider a time-harmonic plane wave 

incident on the periodic inhomogeneous layer from the top region := {x G M'^ | X3 > 6}, where 
d = (ai, Q2, — /3) = (cos 9i cos ^2, cos 9i sin 82, — sin 9i) is the incident wave vector whose direction 
is specified by 9i and 62 with < < tt, < 02 < 2tt and the vectors p and q are polarization 
directions satisfying that p = y^jL/e{q x d) and qJ-d. The problem of scattering of time-harmonic 
electromagnetic waves in this model leads to the following problem: 

curl curl - A;^^ = in X3 > b, (2.2) 

cmlcmlE - k'^qE = in Qb, (2.3) 

u X E = on To, (2.4) 

E = E' + E' in Mi^, (2.5) 

where k = ^/eojloj is the wave number, q(x) = ^(e(x) + i-^) is the refractive index and u is 
the unit normal at the boundary. 

Set a = (ai,a2)0) G and n = (ni,n2) G Z^. The periodicity of the medium motivates 
us to look for a-quasi-periodic solutions in the sense that £'(xi, X2, X3)e~*"'^ is 27r periodic with 
respect to xi and X2, respectively. Since the domain is unbounded in the X3-direction, a radiation 
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condition must be imposed. It is required physically that the diffracted fields remain bounded 
as X3 tends to +00, which leads to the so-called outgoing wave condition in the form of 

E'{x) = ^„e^(""-^+^"^^\ X3 > b, (2.6) 

where a„ = (ai + rii, 02 + ?t-2) 0) G M^, En = {En \ En \ En^) G are constant vectors and 

/3 = / ~ l«riP)5^ if |an| < k, 
[ i(|a„p - if |an| > k, 

with i"^ = —1. Furthermore, we assume that /3„ 7^ for all n G Z^. The series expansion in (|2.6|) 
is considered as the Rayleigh series of the scattered field and the condition is called the Rayleigh 
expansion radiation condition. The coefficients En in (j2.6p are also called the Rayleigh sequence. 
From the fact that div E*^ (x) = it is clear that 

an ■ En + PnE^^^ = 0. 

The direct problem (DP) is to compute the total field E in ]Rij_, given the incident wave 
the refractive index q{x) and the boundary condition on Fq. Since only a finite number of 
terms in (j2.6p are upward propagating plane waves and the rest is evanescent modes that decay 
exponentially with distance away from the periodic medium, we use the near field data rather 
than the far field data to reconstruct the refractive index q{x). Thus, our inverse problem (IP) is 
to determine the periodic medium q{x) from a knowledge of the incident wave E'^ and the total 
tangential electric field x on a plane F^ = {x G M'^ | X3 = a} (a > 6) above the layer. 

2.2 Quasi-periodic function spaces 

In this section we introduce some function spaces needed for the scattering problem (|2.2p - (j2.5p . 
These spaces will play a crucial role not only in the study of the direct problem but also in 
the inverse problem. In [H [121 123) the authors always seek the //^-variational solution for 
the magnetic field based on the facts that the magnetic permeability > is a constant 
and that any vector field H G L'^{Df satisfying that V x G L'^{Df and V • F G L'^{Df 
belongs to Hi^^{D)^ for any bounded domain D C M'^. In this paper, based on the classic Hodge 
decomposition, we are interested in weak solutions in //(curl) of the problem (|2.2p - (l2.5p . that 
is, both E and x E belong to Lf^^(M.^Y. This allows us to solve the scattering problem in a 
general case when ^ is a periodic variable function other than a constant. 

The scattering problem can be reduced to a single periodic cell. To this end, we reformulate 
the following notations. 

Tfe = {3:3 = 6 1 < xi, X2 < 27r}, Q;, = {x G | X3 < 6, < xi, X2 < 2tt}. 

We also need the following scalar quasi-periodic Sobolev space: 

H^inb) = {n(x) = Unixs) exp{ian • x) | tt G L2(0fe), Vn G (/^^(J^b))^ G C}. 

Denote by H^{Tb) the trace space of H^{Q,b) on Ff, with the norm 

\\Khr^=Y. l/n|'(l + |anP)i f^H-Hn), 
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where /„ = (/, exp(ia„ • x))^2(p^) and write H 2 (F^) = (//2(rft))', the dual space to H^iT^). 
We now introduce some vector spaces. Let 

i?(curl , Ob) = {E{x) = ^ En{x2,) exp(iQ„ • 2;) | E„ G C^ 

with the norm 

ll^llH(curl,Q,) = ll^lli2(o,) + l|curlE||i.(^^). 

For x' = (xi, 2:2, 6) G Ffc, s G M define 

i?/(Ffe) = {E{x') = EneMian • x') | K G C^ 63 • ^ = 0, 

nGZ2 

neZ2 

iJ/(div,Ffe) = {E{x') = EneMian • x) | G C^ 63 • = 0, 

neZ2 

ll^lln^div.r,) = E (1 + + ■ a„P) < +00} 

neZ2 

i?ncurl,F6) = {£;(x') = EneMian • x') | £ C^ 63 • ^ = 0, 

neZ2 

neZ2 

and write (Ffe) = i?°(rfe). Recall that 

H;^/\dw , n) = {63 X ^Ir, I E G //(curl , 0,,)} 

— 1/2 

and that the trace mapping from //(curl , 0^) to (div , F;,) is continuous and surjective (see 

[8] and the references there). 

It is well-known (see [20|) that the free space a-quasi-periodic Green function for the 
Helmholtz equation (A + k'^)u = in is given by 

G{x,y) = -^Y T^exp(ia„ • (x - y) + i/3„|x3 - ysD- (2.7) 

^ neZ2 ^^"^ 

We assume throughout this paper that q satisfies the following conditions: 
(Al) q G C^(Ofe) and q{x) = 1 when X3 > b; 

(A2) Im [q{x)] > for all x G Ofc and Im [(?(xo)] > for some xq G 0;,; 
(A3) Re [^(x)] > 7 for all x G Of, for some positive constant 7. 
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3 A quasi-periodic boundary value problem 



Before studying the original problem (j2.2p - (j2.6p . we first consider the following quasi-periodic 
boundary value problem in 0;,: 

cuilcmlE - k^q{x)E = in (3.1) 

ux E = on To, (3.2) 

uxE = f on Tfe, (3.3) 

where / G H^^^T,). 

Lemma 3.1. // the conditions (Al) — {A3) are satisfied, then there exists a unique solution 
E G H{cutI , Qb) to the problem (TO]) - (fg3|l such that 

div \ 

where C is a positive constant independent of f . 

Proof. We first prove the uniqueness part. Let / = 0. Multiplying both sides of (13. ip by E it fol- 
lows from Green's vector formula, the quasi-periodic property of E and the boundary conditions 
(1321) and (1331) that 



/ [|curl^|2 - A;2g|^|2]da; = 0. (3.4) 
Take the imaginary part of the above equation and use the assumption on q{x) to find that 

/ \E{x)\^dx = 0, 
Jb,{xo) 

where B^{xq) C 0;, is a small ball centered at xq with radius e. Thus E{x) = in B^{xq). By 
[9l Theorem 6] we have E G (//^(il^))^. Thus, by the unique continuation principle (see |2H 
Theorem 2.3]) we have = in ilj. 

We are now in a position to prove the existence of solutions. For any V G if (curl , fif,) such 
that X ii^ = on Lq U T^, multiplying both sides of (j3.ip by V yields 

/ [curl curl F-/i:\-E;-F](ix = 0. (3.5) 

There exists at least one element W G -ff(curl , O^) satisfying that vxW = ^ owTq and vxW = f 
on Lft. Then the equation (|3.5p can be rewritten as 

/ [curl [E-W)- curlF - k'^q{E - W) ■ V]dx = - [ [curl • curlF - k'^qW ■ V]dx. 

Let X :={U e iJ(curl , n^), x [/ = on Lq U rj. Then U := E-W € X. Thus the problem 
(j3.ip - (j3.3p is equivalent to the following variational problem: Find U £ X such that for any 

Vex, 

[ [curlU -cuilV -k'^qU ■V]dx = FwiV), (3.6) 
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where F]^[V) = — [cmlW ■ curl V — k'^qW -Vjdx. The proof is broken down into the following 



steps. 

Step 1. To establish the Hodge decomposition: 



X = Xo® VS, (3.7) 
where S = {p £ H^{nb), p = on Tq U Tb} and Xq = e X \ [ q{x)^ • = 0, Vp E S}. 



For U,V eX define 



a{U, V)= [ [curl U ■ curl V - k^qU ■ V]dx. 
It follows from the assumptions (A1)-(A3) on q{x) that 

\a{Vp,Vp)\ >k' [ Re[q{x)]\Vp\'dx > k^7\\Vp\\l,^^^^ = k''j\\Vp\\j,^^^^, 



Thus, for every E G X there exits a unique p £ S such that a(Vp, Vg) = a{E, Vg) for all q €z S. 
Let ^ := E — Vp. Then it is easy to show that E, G and Xq n = 0, which implies the Hodge 
decomposition (13. 7p . 

Step 2. To prove the existence of a unique solution U £ X to the problem f (j3.6p . 
By (j3.7p we may assume that U = ^ + Vp, V = r] + Vq with G and p,q £ S. Then 
the problem (|3.6p becomes the following one: Find ^ G Xq and p £ S such that 



a(Vp, V(7) + a(^, r?) = Fw{Vq) + Fwiv)- 

Since a(-, •) is coercive on VS, there exists a unique p £ S such that 

a(Vp,Vg) = FvK(Vg) VgG5 

with the estimate | |Vp| | /^(cnr/,r2i,) ^ C| I (cur/, ri;,)- It remains to find ^ £ Xq such that 
"^(C)^) = -^iy(^) for a-ll f] G ^0- The bilinear form a(-, •) can be decomposed into the sum of the 
following two forms: 



curl ^ • curl r] + • rjdx, 



«2(e,^) = -k^ I {l + q)^-fjdx. 

JQ.b 

Obviously, ai(-,-) is coercive on Xq, and it follows from [3, Lemma 3.2] that Xq is compactly 
imbedded into Thus, by the standard Fredholm alternative theory there exists a 

unique ^ G Xq satisfying that a{£,,r]) = Fy[r[r]) for all r] £ Xq. Furthermore, | |C| l/i'(c«r/,nb) — 

C\\W\\H(curl,nt)- 

Step 3. To establish the estimate p.4p . 

By Steps 1 and 2 we know that E = + Vp + G //(curl , Qf,) is a solution to the problem 
(I3.1D-(I3.3D with the estimate 



\E\\H{curl,ni) < IIClU(curl,C6) + W"^ P\\ H (curl ,nb) + 1 1 ^^1 U(curl ,^6) < C*! I W^l U(curl .H;,) • (3-8) 
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Recalling that 



ll/IL-i/2/r ) = inf{ll^llH(curi,f76) \ iy X W = onFo and ly X W = f on Fb}, 

it follows from ^ that ||^|U(curi,n,) < C||/||^.i/2(r,). □ 

For / € H^y^{Ti,) define the operator T by 

r(/) = X (curl£; X z^) on Tfc, 

where E solves the quasi-periodic boundary value problem ()3.ip - p.3p . By Lamma I3.H the 
operator T is well-defined. Note that T{f) belongs to the dual space (H^^^"^ (Th))' = -f^curi'^(^^) 

— 1/2 

of ^(jjy (r^) with the duality defined by 

<T{f),g>= / [cutI E -cmlV -k^qE -Vjdx 

— 1/2 

for g G -fffjiv i^b), where V E -ff (curl , ilf, ) satisfies that uxV = gonPf, and uxV = on Fq . The 
operator T can be considered as a Dirichlet-to-Neumann map associated with the problem (13. 1^ - 
(j3.3p and depending on the index q{x). Under the assumptions (A1)-(A3), the above definition of 
r(/) is independent of the choice of V and therefore T : H^^.l^^iVb) {H'l^^iTb))' = H'^/^'^iVb) 
is well-defined. Moreover, it follows from the above equality and Lemma |3. II that 

ll^(/)llji^-i/2(rj,) ^ C\\E\\H{cnri,nh) < ^"■^"^-^''^{rb)- 

This implies that T is continuous from H^^^'^iFi,) to H^^^^"^ {Ff,). 

4 Solvability of the scattering problem 

In this section we will establish the solvability of the scattering problem ()2.2p - ()2.6p . employing the 
variational method. To this end, we propose a variational formulation of the scattering problem 
in a truncated domain by introducing a transparent boundary condition on F^. The existence 
and uniqueness of solutions to the problem will then be proved using the Hodge decomposition 
together with the Fredholm alternative. 

4.1 Transparent boundary condition and variational formulation 

Let x' = (xi, X2, 6) € Ffo for 6 > 0. For E £ (div , F^) with E{x') = Yln&z'^ exp(zQ„ • x'), 
define U : //^^ (div , F^) ^ H^^cml , F^) by 

{nE){x') = (63 X curl£;) X 63 on Ff,, (4.1) 

where E satisfying the Rayleigh expansion condition (|2.6p is the unique quasi-periodic solution 
to the problem 

curl curl ii^ — k'^E = for X3 > b, u x E = E{x) on F;,. 
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The map TZ is well-defined and can be used to replace the radiation condition (j2.6p on Th- Then 
the direct scattering problem (I2.2p - (l2.6p can be transformed into the following boundary value 
problem in the truncated domain Qi,: 



curl curl E — k qE 
V xE 

{cml E)t -n{e^ X E) 



in fift, 
onFo, 

(curl E')t - 7^(e3 x E') on F;,, 



(4.2) 
(4.3) 
(4.4) 



where, for any vector function U , Ut = {y x U) x v denotes its tangential component on a 
surface. The variational formulation for the problem (j4.2p - (j4.4p can be given as follows: find 
E ■.= {E e iJ(curl , Ofe) I 1/ X S = on Fo} such that 

B{E, if) := / [curls • curl^ - k^qE • ^] - / 7^(e3 x E) ■ (es x (f)ds (4.5) 

= [ [{cml E')t - 7^(e3 x E')] ■ (es x ^)ds 

for all ip £ X. We have the following properties of TZ: 

1 _ 1 

1) TZ : (div , F^) (curl , F^) is continuous and can be explicitly represented as 



(nE)(x') 



E 



ne2 



i(3„ 



{an ■ En)an 



exp(za„ • x'). 



(4.6) 



2) Let P = {n = (ni, 712) G | /3„ is a real number}. Then 

1 



Re < TZE, E > 



4.^ J2 



l/3n| L' 



k'^\En\'^ -\an- EJ^ 



-Re<TZE,E> > Ci||divS||^^^_i/2 - C2II-EI 



where Ci and C2 are positive constants and < •, • > denotes the inner product of L^{T 



(4.7) 
(4.8) 



3) 



Im < TZE,E >= 47r2 ^ 



neP 



k^lEr, 



> 0. 



(4.9) 



The representation (j4.6p of TZ can be computed directly from its definition (I4.ip (see [T]) and 
the properties (|4.7p - (|4.9p can be easily obtained using this representation. Furthermore, there 
exists a C > such that for every r] > and E S //(curl , Q^) we have (see [1]) 



Let 



Wu X i^||^-i/2(r^) < C [7?||curli?|U2(n,) + (1 + 1/??)||S| 1^2(0,)] • (4.10) 

S = {p£ HHna)\p = OonTo} 
Xq = {E £ X \ B{E,Vp) = Oyp£ S} 

Then in a completely similar manner as in the proof of Lemma 14. 2| we can establish the Hodge 
decomposition X = Xq © \/S. 



9 



4.2 Solvability of the direct scattering problem 

Lemma 4.1. The bilinear form B(-, •) defined by ( |^.5[ ) is strongly elliptic on Xq, that is, for all 
Wo G ^0, 

ReB{wo,wo) > C\\wo\\x - p{wo,wo) 
for some constant C > and a compact bilinear form p{-, ■). 
Proof. Let M be a positive constant to be determined later and let 

6i(u;o, 9?o) = / [curltDo • curl^ + (M - k^q)wQ ■Tp^'\dx - / 7l(e3 x wq) ■ (es x ■^o)^'^ 
b2{wo,(fo) = -Ml wo-Tp^dx. 



Then B{wo, (fo) = bi{wo, ipo) + 62(^)0, (fo) for wq, ipo £ Xq. By the properties of TZ it follows that 

—Re < 7^(63 X Wo), es xwo > 

> Ci||div(e3 X u;o)|l5^-i/2(p^) -C2||e3 x t«o||^^-i/2(p^) 

> Ci||div(e3 X u;o)||^-i/2^p^^ - C3^^||curlu;o|li2(n,) - ^3(1 + ;^)^lko|li2(n,), 
where Ci, C2 and C3 are three positive constants and 77 > is arbitrary. Thus we have 

Rebi{wo,wo) > ||curlu;o||i2(Q^) + (M- /c2goo)||u'o||i2(f^^) 

-C3r/2||curlu;o||i2(n,) - Csil + l/r]f\\wo\\l2^^^^ 
= (1 - C37?2)||curl«;o|li2(j,^) + [M - fc^^oo - ^3(1 + l/vf]\\wo\\l2^^^-,, 

where goo = ^^^xeR^ < Choose ry sufficiently small and M sufficiently large so that 

Rebi{wo,wo) > Co(||curl?i;o|li2(n,) + lko|li2(f^^)) (4.11) 

for some constant Cq > 0. This, together with the fact that Xq is compactly imbedded in 
(^■^(f];,))^, yields the desired result. □ 

Theorem 4.2. Assume that the conditions {Al) — {A3) are satisfied. Then the problem l\2.^ — 
l\2.6^ has a unique solution E G i^joc (curl , such that 

\\^\\Hioc{cMr\,M.\) ■= max||£'||^(cui.l,Qa) ^ C'll^'llH(curl,Q(,), 

where C is a positive constant depending on the domain and q. 

Proof. It follows from Lemma [4.1l and the proof of Lemma [3. II that there exists a unique solution 
E S -ff(curl , Qft) satisfying that | l-El |/f(curl ,^2;,) ^ C||-^*llH(curl,Oi,)- h remains to extend E{x) to 
be a function in Hioc{curl ,R\). Suppose 63 x {E - E^)\r^ = EneNxN ^ne'""'^ e if-^/2(div , Fb). 
Let 

E'{x)= (^„xe3 + i?„e3)e*""-"+'''"("^-'), x^ > b 

nSNxN 
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and let E'^ satisfy that div E^{x) = for X3 > b. Then we have i?„ = ^(es x An) ■ an- Thus 



nSNxf 



AnX 63 + ^(es X An) ■ a„e3 

Pn 



Ja„-x+il3n{x3-b) 



X3 > b. 



Define E{x) = E\x) + E^{x) for X3 > b. Then it is easy to prove that E e if (curl , ^a\^b) with 



\E\ 



< C| I i?* 1 1 //(curl, Hi,) for any a > b, so E ^ H{cutI ,Qa) for any a > b, that 



S e ifioc (curl , M;^) with the required estimate (I4.12p . The proof is thus completed. 



IS, 

□ 



5 The inverse problem 

Let a > b and assume that there are two refractive index functions qi [i = 1,2) satisfying the 
assumptions (A1)-(A3). For g G LliTa) let the incident waves be of the form: 

£"(x,5) = curUcurU / G{x,y)g{y)ds{y), x < a. (5.1) 

Write the scattered electric field and the total electric field as Ef{x,g) and Ei{x,g), respectively, 
indicating their dependance on g and the refractive index function qi (i = 1,2). 

For the refractive index q^ denote by Tj the corresponding Dirichlet-to-Neumann map associ- 
ated with the problem p.ip - (|3.3p with q replaced by (^j (i = 1, 2), as defined at the end of Section 

El 

Lemma 5.1. If Ti{f) = T2U) for all f G i?7^/^(div , F,,), then 

/ Ei{x) ■E2ix)[qi{x) - q2{x)]dx = 0, 

where Ei, E2 £ H{curl,i^h) solve the problem — (jc/.^ with q replaced by qi andq2, respec- 
tively. 

Proof. Let Ei and F2 G //(curl , Qb) be the solution of the problems 

curl curl £^1 — k^qiEi =0 in Vt^, u x Ei = on Fq 

and 

curl curl i<2 — ^^(72-^2 = in 0;,, u x F2 = on Fq, u x F2 = v x Ei on F^, 
respectively. Let E = F2 — Ei. Then it is easy to see that 

curl curl E' — A;^g'2-£' = k'^{q2 — qi)Ei in Ofe, 
1^ X E = on Fo U Ffe, 
v X curl£^ = on F^, 
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where the last quahty is obtained from the assumption Ti =T2- Thus, it follows from the Green 
vector formula that 

/ {q2 - qi)Ei ■ E2dx = I {cuilcmlE - k'^q2E) ■ E2dx 

= To f (cuilE ■ cm\E2 - k^q2E ■E2)dx 

= I (E ■ cuT\cmlE2- k^q2E ■E2)dx 

= ^ [ {E- k'^q2E2 - k'^q2E ■ E2)dx = 0. 

The proof is thus completed. □ 
For g £ L'f{Ta) appearing in the incident waves (jS.ip . we define an operator F : LfiTa) 

Fig) = e3 X E{x,g) on Tfe, 

where E{x,g) solves the problem (|2.2p - (j2.5p with the incident wave E^{x,g). The operator F 
can be considered as an input-output operator mapping the sum of the electric dipoles to the 
tangential component of the corresponding total field on Ffe. Moreover, for all g £ Lf{Ta), the 

— 1/2 

operator F has a dense range in (div , Ffe), as stated in the following lemma. 

— 1/2 

Lemma 5.2. The operator F has a dense range in (div,Ff,). 

Proof. We only need to prove that F* : ^^^(curl,F5) — J- L^(Fa) is injective. First, we show 

— 1/2 

that for any f £ (curljF;,), F*{f) is given by 



curlycurlj^ / G{x,y)cuTl{V~^{x) — W{x)) x e^dsix) 



(5.2) 



where the superscripts + and — indicate the limit obtained from M3\il;, and Qb, respectively, 
and for any a > b the function V E ff(curl , Q^) n ff(curl , Qa\^b) solves the problem 

(5.3) 
(5.4) 
(5.5) 
(5.6) 
(5.7) 

and satisfies the Rayleigh expansion condition (|2.6p with a replaced by —a for 2:3 > b, that is, 

V{x) = y^eiK-^+/3;^3)^ ^3 > 5 (5 8) 

with a'n = (— ai + ni, —02 + ^2, 0) G R^, Vn £ and 



curl curl V — k^V 


= 


for 


X3 > b 


curl curl V — k^qV 


= 


in 




uxV 


= 


on 


To, 


V X y+ — V X - 


= 


on 


n, 


[curl — curl 


= 7 


on 





(/c^ - |a^P)2 if la'^l < k, 
i(|a^P - A;^)2 if \a'^\ > k. 
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In addition, the function W is given by 

W{x) = Ke*«""-^+'^"(26-x3))^ 3.3 < 5. (5 9) 

In fact, for any f £ ^^^(curl , F;,) and g £ "^^^(div , Ff,) we have 

<i"5,/>^-i/2(di,_r,)xH-i/^(curi,r,) 
= ux E{-,g) -Jds 

= u X E{-,g) ■ [curiy+ — cmlV^]ds 

= / [{i^ X E ■ curiy+ - 1/ X y+ • curlS) - {v x E ■ curlF" -uxV~ ■ c\iT\E)]ds, 



where the transmission conditions (j5.6p and (|5.7p have been used. It follows from the Maxwell 
equations ()5.4p and (|2.3p and the boundary conditions (|5.5p and (|2.4p that 

/ [i^ X E- curly- -uxV~ ■ curl E]ds = 0. (5.10) 

On the other hand, from the Rayleigh expansion conditions (j2.6p and (jS.Sp it is derived that 

/ [ly X E ■ curl y+ — u X ■ curl E]ds 

= I [{ux E' ■ curl -uxV^ ■ curl + {u x E' ■ curl 1/+ - x 1/+ • curl 

[vxE' ■ curiy+ - X y+ • curl£;^]ds. (5.11) 



Similarly, from the definition of E^ and the Rayleigh expansion condition (|5.9p it follows that 

/ [i/x£;*-curlW-z^xiy-curl£;*]ds = 0. (5.12) 
The equations (|5.10p - (|5.12p together with the fact that V = W onVi, yield 
<Fg,f > = I [uxE' ■ curiy+ - u x V+ ■ curl£;*]cis 
= [ [u X E' ■ cmlV"^ - u xW ■ curl E']ds 
= I [u X E' ■ cmlV^ - u X E' ■ cmlW]ds 
= j ux E' ■ (curl y+ - curl W)ds. 
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Substituting the expression (|5.ip of E"^ into the above equation and exchanging the order of 
integration we get 



< Fg,f >= / g{y) • curlj^curlj^ / G{x,y)cuil[V^ (x) - W{x)] x esds{x) 
which imphes ([57 



ds{y), 



— 1/2 

We now prove that is injective. Suppose F*(f) = for some f e H^ ''"(curl,r6). Define 
Uhy 



U{y) := curlj^curl. 



G{x, y)h{x)ds{x) 



y G M3\rfe, 



where h = curl (y+ — W) x 63. Then 63 x U{y) =0 on Fq. It is clear that U{y) is a —a- 
quasi-periodic function satisfying the Rayleigh expansion condition (j2.6p when 2/3 > a. By the 
uniqueness of solutions to the exterior Dirichlet problem (see [2j) we have U{y) = when 7/3 > a, 
which together with the unique continuation principle ([llj) implies that U{y) = when y^ > b. 
Now from the jump relation 63 x U^{y) — 63 x U^{y) = on and again the uniqueness of 
solutions for the exterior Dirichlet problem for y^ < b we get that U{y) = when 1/3 < b. Thus, 
h(y) = 63 X curl [[/"'" (y) — U~{y)] = on F^, which, together with (jS.Sp and (j5.9p . implies that 

63 X y+ = 63 X VF, 63 X curiy+ = 63 X curlVF on r^. (5.13) 

Since V and W satisfy the Maxwell equation curl curl — k'^E = in the regions X3 > 6 and 
x-^ < b, respectively, then it follows easily from the transmission condition (|5.13p and the Rayleigh 
expansion conditions (15. Sp and (15. 9p that V = for X3 > b and W = for X3 < b. Thus, by 
(j5.6p we have x = on F^, so F G i7(curl , fif,) satisfies the problem (j3.ip - (|3.3p with / = 0. 
By Lemma l3.ll we have V = in ri^. Thus, / = [curly — curll/ ]t = 0, which completes the 
proof of Lemma 15.21 □ 

Combining Lemmas 15. II and 15. 2[ we have the following orthogonality relation for two different 
functions qi (i = 1,2). 

Lemma 5.3. Let the incident waves E^{x,g) be defined by (j5.jp . // 

63 X Ei{x,g) = 63 X E2{x,g) on F^ (5.14) 
for all g G L|(Fa) and some a > b, then the following orthogonality relation holds: 

Ei{x) ■ E2{x){qi{x) - q2{x))dx = 0, 



where Ei, E2 & H(cml,^}b) solve the problem i\3.1\i — (jg.gp with q replaced by qi and q2, respec- 
tively. 

Proof. From the equation (|5.14p . the uniqueness of solutions for the exterior Dirichlet problem 
and the unique continuation principle it follows that Ei(x,g) = E2{x,g) for all x-^ > b. This 
implies that 

63 X curl£'j^(x,5f) = 63 X curl£'^(x,5) on F;,. 
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Since [es x curl Ej~ {x; g)]\r^ = for j = 1,2, then we have 

63 X curlE^{x,g) = 63 x curl£'^(x,5) on Th. 
By the above two equahties and the definition of Tj we have 

Tiies X Ei{x,g)) = T2{e^ x E2{x,g)) 
for all g € L^(ra). The continuity of Tj (j = 1, 2) and Lemma [5.21 lead to 

n{f) = T2{f) yfGH;'^\diy,n). 

This together with Lemma l5 . 1 1 gives the desired result. □ 

We are now ready to prove our main result for the inverse scattering problem. 

Theorem 5.4. Let qj {j = 1,2) satisfy the assumptions (Al) — {A3) and let qj depend on only 
one direction xi or X2 with j = 1,2. If 

63 X Ei{x,g) = 63 X E2{x,g) on Ta 

for all g G LliTa) with some a > b, where Ej{x,g) solves the problem {\2.^ — {\2.5\i with q = qj 
(j = 1,2) corresponding to the incident wave E^{x,g) given by ()5. ip . then qi = 52- 

Proof. By Lemma [5.31 we have the orthogonality relation: 

/ Ei{x)-E2{x)[qi{x)-q2ix)]dx = 0, (5.15) 

where Ei, E2 G //(curl , ilj,) solve the problem (|3.1|) - (|3.3|) with q replaced by qi and q2, respec- 
tively. 

We now look for solutions to the problem (I3.ip - (l3.3p in the following form: 



E{x) = (0,0,^3(xi,X2)) = {0,0,v{xi)u{x2)) 

with the scalar functions v and u satisfying the following quasi-periodic conditions: 

f (xi)e2^°i^ = v{xi + 27r), u{x2)e^'''^'' = v{x2 + 2^). 

It is clear that such a function E is a-quasi-periodic and satisfies the boundary condition (|3.2p . 
Without loss of generality, we may assume that qj{x) = qj{xi), that is, qj depends only the 
j;i-direction with j = 1, 2. Substituting such E into the Maxwell equation (j3.ip and noting that 
curl curl = —A + V(V-), we find that 

v"{xi)u{x2) + v{xi)u"{x2) + k'^q{xi)v{xi)u{x2) = 0, xi, X2 G (0, 27r), 

which implies that 

+ k^q{xi)v{xi) = — ^ = A 
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for some constant A, where xi,X2 G (0,27r). Following the idea of Kirsch [17] . we construct a 
special kind of solutions v by considering the following quasi-periodic Sturm-Liouville eigenvalue 
problem: 

{v"{xi) + k'^q{xi)v{xi) = Xv{xi), xi £ (0, 2tt) 
w(xi)e2*"i'^ = v{xi + 2it), 
v'(xi)e2*"i^ = t;'(xi + 27r). 

The eigenvalues A^, and the corresponding eigenfunctions Vn, normalized to Vn{0) = 1, have the 
following asymptotic behaviors as n — t- oo (see |29|): 

2 f2TT 



v^{xi) = exp 



i{±n + —)xi 
zvr 



+ 



which are uniform in xi G [0, 27r]. We also consider the following quasi-periodic boundary 
problem for u: 

(m ■ / ""(^2) - Ku{x2) =0, X2 e (0,27r) 

^ ' ' \ u{x2)e^'°"^'' = v{x2 + 2-k). 

The non-trivial solutions to the problem (II) can be written explicitly as 

Un{x2) = c„,ie^""2 + c„,ie-^""2 , A„ / 0, 
where Cn,i and c„^2 are constants satisfying 

Cn,l = Cn,2 (e^^-v^n - e^^-"^) / (e*2™2 _ g2^V^„\ (5 ig) 



Now, let E^^ = v^{xi)u^{x2) be the third component of = (0,0, i?^^) corresponding to 
qi{xi) and let E^^ = v^{xi)u^{x2) be the third component of E^ corresponding to q2{xi). It 
follows from (jS.lSp that 

0=/ E^^n{xi,X2) ■ E^^rn{xi,X2)[qi{xi) - q2{xi)\dx = bAl''^A2''"^ , (5.17) 



b 

where 

f27r 



^i'"': = ^"fe(xi)-g2(^i)]e^("-"^)^^dxi + o(i)+o(l), 

and Crtj, Cm J satisfy (j5.16p with j = 1, 2. For arbitrarily fixed / G N, letting m = n — I gives 
^r''" = j\q,{x,)-q2{x^)y'^-d{x,)+o{£^, 
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We can always choose appropriate constants Cm,2 and Cm,i satisfying (j5.16p such that ^^+''™- ^ q 
for sufficiently large m. In fact, we may assume that I is a positive number since other- 
wise we can take n = m — I' for some positive V instead of /. Now choose Cm,2 = e^'^^^^ . 
Then, by (|5.16p . \cm,i\ > Ci for large m with some positive constant Ci independent of m 



and 



^2 I 

m —7- +00. Letting m — )• +00 we conclude from (|5.17p and the above discussion that 

''"(91(2:1) -g2(xi))e^'"idxi =0 
/o 

for every / G N, which implies that qi = q2- The proof is thus completed. □ 



Jq^ Cm,2e '^'^^'^"^^^dx2 tends to +00 as m — )• 00. This implies that \Al^ — )• +00 as 
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